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ABSTRACT 


PitemaschuccuTe Of barotropically unstable disturbances 
in the tropics is studied with a two-level quasi-geostrophic 
model. An Ekman layer is attached to the lower boundary. 
The equations are linearized and the most unstable mode is 
found numerically by use of the initial value technique. 
CISK heating, horizontal friction and shear are introduced 
to give more realistic solutions. Computations are made for 
Pwo shear zone wind profiles. The first is represented by 
U, = U2 = -U,tanh y/y, and the second is a mirror reflection 


of the first and is given by Us = The model is run 


iL 
for 142 days after which growth rates are determined as a 
function of the non-dimensional wave number KYo° The 
associated phase angles and amplitudes for the wave numbers 
with maximum growth rates are given. The solutions with 


heating do not have a strongly preferred mode of maximum 


growth. 
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i. INTRODUCTION 


synoptic scale disturbances in the tropical atmosphere 
have been recognized for a long time (Palmer 1951, Riehl 
1954). This paper is concerned with wave disturbances that 
occur in the lower troposphere and makes no distinction 
between easterly waves and equatorial waves. Starr and 
Wallace (1964) showed that in the lower troposphere these 
disturbances have a "cold core" structure; the rising air is 
cooler than the air which is subsiding. Warm core structures 
at higher levels with significant correlations between 
memperature and vertical velocity have been found by Chang 
et al. (1970) and Nitta (1970). Vertical velocities associ- 
ated with these disturbances in general are rising air to 
the east of the wave axis and subsidence to the west (Yanai 
ema Nitta, 1967). 

The energy source for these waves is not known. 
Bemeoclinic instability is unlikely because of the smallness 
of the Coriolis parameter. Charney (1963) performed a scale 
analysis which indicated that, in the absence of condensation, 
large-scale motions in the tropics are governed by the baro- 
tropic vorticity equation. Charney (1963) suggested that 
easterly waves arise from barotropic instability. The 
spectral analysis of Wallace and Chang (1969) indicate a 
conversion of kinetic energy from the mean flow to the 


disturbance. 





A third possible energy source for these disturbances 
Homlaeenu Neat of condensation. 

This mechanism is consistent with observations of warm 
core structures at upper levels (Chang et al. 1970). This 
heating has been parameterized by Ooyama (1964), Charney 
and Eliassen (1964), and Kuo (1965) in such a manner that 
moe condensation heating is proportional to the convergence 
of moisture in the Ekman layer. This leads to a simple 
i@avabllity which is called Conditional Instability of the 
Second Kind (CISK). Yamasaki (1969) studied the influence 
of CISK on linearized disturbances in an atmosphere with 
vertical shear. The solutions displayed many features that 
are observed in actual tropical wave disturbances. Holton 
(1971) has obtained similar results with a model that con- 
mens a2 mOving heat source and vertical wind shear. These 
studies strongly suggest that the tropical disturbances are 
maemver by condensation heating. It is known that the CISK 
Meenanism must have a pre-existing disturbance to organize 
Bie convective elements before CISK can Operates has. bLe— 
existing disturbance could arise from barotropic instability 
as was demonstrated by Babee (1970). Kuo (1949) examined the 
eeao1 lity characteristics of a barotropic zonal current and 
momma that instability is possible only when the absolute 
vorticity gradient reverses sign in the region. 

The purpose of this thesis is to extend the work done 
Bema liijiams et. al. (1971) to include cases with CISK heating, 


horizontal friction and vertical shear. Since the Williams 





et al. (1971) paper is referenced many times, it will 


hereafter be referred to as WSN. 





JI. THE FORECAST EQUATIONS 


The quasi-geostrophic prediction equations will be used 
as the forecast equations following WSN. It is recognized 
that the quasi-geostrophic approximation is not very accurate 
near the Equator, but Matsuno (1966) has shown with a baro- 
tropic model that quasi-geostrophic motions are identifiable 
mearm the Equator if the wavelength is not too long. 

Yamasaki (1969) has found that the error is small at 20 
Peam@@ees latitude in a study with a baroclinic model. 

Phe wdUasi—feosurophic equations are simplified through 
Pmenincroduction of a two-level model which was also used 
by Bates (1970). The equations are linearized and solved 
mummers cally by the initial value technique. The growth rates 
and wave structures are then computed. 

A simple two-level model is constructed by dividing the 
memes phere invo four layers with a constant pressure differ- 
ential of Ap/2 (Fig. 1), numbered 0 to 4 from top to bottom. 
Mire DOundary condition at p = Dy is WO, = Oreeworneeho ae rOh- ors 
vertical propagation of energy. Charney (1969) has shown 
that such energy propagation is unlikely. At p= p, the 
boundary condition is w = We where W Bom ne wire bLOlal iy 
mieiiced @ at the top of the Ekman layer. Charney and 
Eliassen (1949) derived the following expression for w 


S 


A 
O, = “Py (se)? sin Cat ; C2) 
fe) 





where A. is the Kinematic eddy viscosity, ais the surface 


inflow angle, and oH 1s the surface geostrophic vorticity. 


p Level 
Py W=0 0 
¥y 
1 
Ww 
Py zeae hy 6. - P, 
Ve 
3 
W= Ww 
Po ~ y 


Fig. 1 Two-level prediction model 


It has been shown by Holton et al. (1971) and Chang (1973) 
that this expression is reasonably valid as far equatorward 
as 10 degrees, but may be invalid when close to the Equator 
Memsoutcm Of the ITCZ. The surface geostrophic vorticity can 


be approximated by 


2 
Cy = $3 oe We ¥2 ° (202): 
The quasi-geostrophic vorticity equation in pressure 


coordinates with horizontal eddy diffusion is 


ow HT 


S 9p + kxp-V(V) + Bo = ~f¢ Beawiy. (2.3) 


oe O dp 


Observe that when w = 0 and A, = 0, (2.3) becomes the baro- 
imeontc VYOrticity equation which is sufficient to describe 


the main instability. 


1G, 





Apply this equation at levels 1 and 3 which gives 


aw W) 

,) 2 2 1 ae 4 

ewe a Bo ox > fo Bp Aa’ Yi > (82-7) 
ow (w)—-w, ) 

oe ae 2 3 eee 4 


If it is desired to stop at some intermediate point, say the 
MmeemOopailsce., ravher than to include the Entire atmosphere, 
smaller layers could be used. 

Wext consider the quasi-geostrophic first law of thermo- 


dynamics in the form 


2 
-uU KQ 
o.060OoOW ow 0 2 
ae ae + kKxVP (=) + — w = ——= (2.6) 
aul-op dp Ps of. 
where 
_ _ | 
-~RT , wf, & 
Oma" = [ Re die (2.6a) 
P& p 


and T is the temperature obtained from the standard atmosphere. 
The heating follows the formulation of Ooyama (1964), 
Charney and Eliassen (1964), and Kuo (1965), i.e., the 
condensation heating is proportional to the convergence of 
moisture in the Ekman layer. The form used was given by 


Ogura (1964) 


af 





Cad 


ais omer Co) 
Q oe ees Vere 
Ee 2 ry 


Nh Ww 9 
dp 4 


where n is a non-dimensional parameter. This relationship 


(25 


mecca 1n both rising and subsiding air since a Fourier 


@eempanslon in x is used; therefore the amplitude of heating 


is divided by two. 


When (2.6) is applied at level 2 the result is 


2 
Uy Ns APO 505 = KQ> 





3 - 
ge V-¥3) + kxV(—3 >) VC 03) - RO FOF 


O O 


Wetane the following quantities 


V1 + v3 
Yu = a on 9 
Ui = Wve 
Yn = aC 
which implies that’ 
¥, = Wy + Ups 
¥2 = Vy 7 Yn 


(2, 


C2: 


(2 


(2 


(2x 


7) 


8) 


9) 


10) 


neil) 


1) 


Here Vn Hee proporcional to the layer thickness and is there- 


fore a measure of the mean temperature in the layer. 


Use 


m@i@ese definitions, ada (2.4) and (2.5) and divide the result 


by two, which gives 


eZ 





oe 62 2 2 
ae Vy + KxVby VOU) + kxVb ns VCVOU,,) 


dW ae) 
M 4 4 


For the second forecast equation, subtract (2.5) from (2.4) 


and eliminate », with the use of (2.8) which gives 


3 (Veuve + XV VOW Ue) Yap + Vy VV Yay) 





O 
dW f_w -u~KQ 
T ej 2 2 h 
° 2 pemey oor. a ry eee (2.14) 
where 
2 
2 21 
un 5 (2.1424) 

p ols 


These are the prediction equations for the model. To 
Hinearize the equations, the flow is separated into an east- 
west current that varies Onyeaaney and aucmall departure 
meom this flow. Thus it is possible to treat waves in x 


independently and the fields can be defined as follows 


Wy = E(y) + ACy,t)cos kx + Bly,t)sin kx, 2 5D 


F(y) + C(y,t)eos kx + D(y,t)sin kx , 225 ALS 


& 
3 
It 


13 





where k is the x wave number. The coefficients A through D 
are Fourier amplitudes of the disturbance and E and F are 
bie coefficients of the mean zonal fields. 

Substitute the expressions for Yur and Vr Om Cee 3) 
eeparave the various sine and cosine terms, and neglect all 
products of the quantities A through D. Equate the coeffi- 


eients of the cosine terms which gives 


3 OCA aD, 59E 9°B , 3F 8°D 3°E 3 3p 
mo TCA) = ame hy ae 
dy y dy dy dy 
2 ,9E oF sc 2 
- k° (x= B+ = D)j - B_Bk - K(—s5 - k”)(A-C) 
oy dy O ay° 
y 2 
oA 2 BoM 4 
“- Ay (—y NE rs ke ae Coma) 
dy dy 
The sine terms give 
9 (8B 2) 2 ype EDA OF OC, PE, RG 
ay aye 2 ay- by? ay? 
+ QE a + 2B 0) + B_Ak - K(oe - x2) (B-D) 
dy dy O ay* 
4 2 
9°B 2 3°B 4 
eek a tc) Ques) 
dy dy 
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Repeat the procedure for (2.14) for the cosine terms which 


gives 
9 ,3°C 2 2 9F 8°D . OF 8°B DE, 2. 2 
ae o> Cin = Cll )) = k [—— Sap POS Ee SS Di 
dy dy oy dy ay oy 
3 3 
Se =) = —_ = aa BJ - 8 Dk 
dy dy dy 
2 y 2 
A) 2 aC 2 Bere y 
+ OS = 5) eer k")(A-C) + A, (— - Be ay Is 3 ear 
dy dy dy 
(2 12) 
The results for the sine terms are 
2 | 
Be _ pce - put) = k Coe ct ee) + OE ack? - 12) 
ot 2 
dy dy dy 
2 2 3 3 
_ 9E aC _ oF oh 5 os A + a C] + 6 Ck 
dy day dy oy dy dy 
2 ' 4 2 4 
Ps) (2 = ke)(B—p) + a.(2 = one MP ak pd) , 
2 H q S 
dy dy dy 
(2220) 
where 
ime ha, 7S : 
ree oo (=) (Oa) 
OR8 e) 
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and 


A 
ia n M,% 1 
S = F Po =) a (2x22) 


mnhese are the same equations used by WSN except that in 
their equations F = 0, n = QO, and AL = 0. 
It is necessary to derive the computational equation for 


vertical motion beginning with 


ow OW dW KQ 
ee ee Ll _t ty ie 
sccm ox’ ’m jy of. 1; (2.23) 
where 
dy 
a aa 
eekatre Sy ; (2.24) 
O 
oy 
See 
ene f, ox (2.25) 


Sieevitute equations (2.24), (2.25) and (2.16) into (2.23), 
collect terms, and again neglect all products of A through D 


which gives 


ne ac ok (3E p _ BF 
Wy = ing {cos kx [5 (sy D ar B )] 
: KQ 
aD k Of of of 2 
+ sin kx [ a ag - “by = xy ADI oe 2F | 
G22 26) 
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Both the disturbance and the average vertical motion are 
Meee oyeunis €duation. The prediction equations and the 
Wercical motion equations are now in the form to be used in 


the model. 


Tey 





IIIT. BOUNDARY CONDITIONS AND FINITE 
DIFFERENCING SCHEME 


The finite difference scheme that was used is illustrated 


below with a sample variable M,: 


aM : 

ay > OH Myay 7 My_y? (3.1) 
9M 1 

- = = (My 44 - eM, F M,_4); G22) 
aM = (mM. .-2M. .-M.)-(M.-2M, .4M. .)] 


(35 3.) 


feere | 1S the grid index and H is the distance between grid 
poncs. 

Centered time differences were used for all quantities 
meeeou Chose involving friction and heating. The friction 
and heating terms were computed at time (t - At). Initially 
a forward time step was used. The second order equations 
for time tendencies were solved by the exact method of 
Heemtmyer (1957, p. 101). 

Yanai and Nitta (1968a) studied the problem of using 
finite difference approximations to solve dynamic instabil- 
ity problems of non-divergent barotropic currents. They 
showed that for either symmetric or anti-symmetric zonal 


G@urrents, the exact boundary conditions can be replaced by 


18 





plactnPp riesid boundary conditions at a distance equal to 
one-half the width of the shearing wind belt. Therefore 
rigid boundaries are placed at y = -W/2 and y = W/2, where 

W is the total width of the computational field. Except 

for the E and F fields, all variables have the boundary 
condition A =Bz=CeH#De=0atye= #W/e. Second derivatives 


miaae are required on the boundary are set equal to zero. 


IY, 





y O 


-4 


“10.0 -5.0 0 


Fig. 2 Mean wind profile U 


20 


2.0 10.0 


ave tanh We 





IV. WIND PROFILE AND INITIAL CONDITIONS 


Woeleanvestvigated Che structure of tropical disturbances 
Meee a Similar model and the wind profile Ui tanh WS 
iemers proposed to expand on that paper to include a sinusoidal 
@m@e@er1 le for an initial disturbance and a thermal wind field 
Which will be a mirror image of the mean wind field used 
by WSN. 

if the waves are to grow, the zonal wind profile must 
be barotropically unstable. Charney (1963) suggested that 
a shear zone could be formed by the transport of angular 
momentum from different regions. This would occur when the 
imeercvropical Convergence Zone is shifted away from the 


fomecor. The unstable profile chosen to represent this is 


TT = |] — =~ | 
Uy U2 U. tanh Wee (4.1) 
which is illustrated in Figure 2. To represent a thermal 


wind field the unstable profile is 


U, = -U. sale tanh Y/Y: Ce) 
Lipps (1970) has used the profile represented by (4.1) 
mma 2nviscid study of barotropic instability in the tropics. 
This wind field is anti-~symmetric with respect to the 


Intertropical Convergence Zone, located at y = 0. 


eal 





For the experiments which have the wind profile (4.1), 
a range of non-dimensional wave numbers (ky) from 0.10 
to 3.00 will be used. This wide range of wave numbers is 
necessary because it is uncertain which wave number will 
be the most unstable mode. Those experiments with the (4.2) 
meee proiile will be run over a range of 0.10 to 1.20. One 
Saeetrament with the latter wind profile containing CISK 
meeoing was run over a wave number range of 0.10 to 3.00. 

Py is introduced into the model as a sinusoidal distur- 
Semce in y and Von Caewemirom Zero as a result of the Ekman 
[Memerton term. CiSK heating Is added which greatly increases 
MEeeeOwoll raves Of all fields, particularly the short waves. 
femeeronoal friction is also added. This represents a mixing 
ereune smaller scale waves and smooths out perturbations 


mie vertical velocity field. Thus at time t = 0 


A = sin = (4.3) 
B = 0 (4.4) 
C= 0 (4.5) 
DS © (4.6) 


E = ioe Uy Wie InLcosh (y/y J or 0 Cave) 


ty 
iH} 


a Ue ee Infcosh Ayia or 0 Cie) 


ye 





V. COMPUTATIONAL PROCEDURE 


The following constants were used in all computations 


Emeeepu those related to heating and friction which were 


zero if the 


At 
n 


90/dp 


heating and/or the friction are neglected: 


iOOmecemua bars. 

400 or 200 or 100 kilometers, 
4000 kilometers, 

LOemeters per second, 

5 x 107° per second, 


O.4 meters per second squared per 
CoMeibar squared , 


0.0 or 2.29 x 10721 


DemmeINeUCe pele sccond. 
SUP eceuuilvars, 

itjemeters perssecond Squared, 

eax 10° Newegemectr second Squared, 

22.5 degrees, 


2.49 x 107° 


peremever squared , 
DUmOnNeouOn Ie.5 kilometers, 
OmOr nas Or O25 hours, 

ZO ES To) «2 allem 0) 


-0.6 degrees per centibar. 


A forecast period of 142 days was used so that the most un- 


stable mode would have reached the maximum growth rate. For 


most cases, it was found that the growth rate at 110 days 


was the same as the growth rate at the longer period. 


oe) 





WSN introduced a random disturbance into the Yu field. 
meewas LOuUnad that by using a sinusoidal profile for the 
Gisturbance, similar growth rates were obtained except in 
ime cases where heating was added. The sinusoidal profile 
was used as the disturbance field for all experiments with 
Variations of U, 8, CISK heating, horizontal and surface 
friction, and the non-dimensional wave number Oey lc Growth 
foeves were then computed by assuming the amplitude could 


be written as 


a, = a, exp(nU./y,) t,; (Soil) 
Qo = a. exp(nU/y,) ty, Cee) 
where n is the growth rate and a, iS a constant. If the 


ratio a5/a, Pomc Lie  gieowln rare can be shown to be 


Veen an7as) 
f= ES Se ao 
Orrec 1 


The dimensional doubling time is 


ie in 2 


T=> ——. Gwe) 
nu. 
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VI. RESULTS 


A series of experiments were conducted over a wide 
meme ty Of conditions. The total distance over the grid 
was 4000 km. A grid mesh size of 50 km and a time step of 
1.0 hours was used for all experiments except two cases. 
Memone Case a2 grid mesh of 25 km with a 0.5 hour time step 
was used and in the other case a grid mesh of 12.5 km and 
a 0.25 hour time step was used. 

The more important experiments are given in Table 1 as 
Pmeemrnction of U, 8, heat, friction, and ky,. The growth 
mewent 15 also given. The table shows whether or not CISK 
Meecing and horizontal friction were included for each 
experiment. All experiments given in Table 1 included 
Semeace iriction. For each experiment, calculations were 
meade to determine the non-dimensional growth rate as a 
function of the non-dimensional wave number KY y S56 thae. ene 
Pesults could be compared ©O those obtained by WSN. The 
Venue of KY, that was associated with the maximum growth 
rate ig shown in the table. 

ieeomould be mentioned that the value used for the non- 
dimensional heating parameter n was too large. A more 
realistic value for n would have been 1.0 instead of 2.0. 
Test runs were made for experiments 1 through 4 with 
nh = 1.0 and in all cases, the growth rates were reduced to 
about one-third the value obtained with the larger value 


of 71). 


25 





Table 1. Compilation of numerical experiments 


eS 


0 
ahi. 


Bee Fe 





Experiment 1 was run with the (4.1) wind profile and 
included CISK heating and surface friction. The beta term 
and horizontal friction were excluded. Garcia (1956) found 
that the shear zone profile (4.1) is barotropically un- 
stable when 0 < KY < 1. Betchov and Criminale (1967) com- 
irea che Gigensolutions for this profile throughout this 
range and found the maximum instability at Le a 0.45 which 
agreed with the results of WSN. The wave number of maximum 
Mepebility (0.40) agrees quite closely with that value but 
the growth rate (n = 0.770) is much greater than the value 


obtained by WSN (n 


0.153) even when the smaller value of 

n was used. Figure 3contains the growth rate profile for 

the full range of wave numbers. The maximum growth rate 
corresponds to a doubling time of 0.42 days and a wavelength 


of 6280 km. The sharp irregularities between 0.10 and 1.00 


26 





pese because the wave number increment in that region is 
wesoeas Compared to 0.50 between 1.00 and 3.00. In the 
non-heating cases of WSN, there was no growth for KY alae 
mien experiment 1 there is large growth which must be due 
Mampare tO the CISK process. The Figure shows that the 
Browth is nearly independent of the scale. Figure 4 contains 


mae amplitudes of Ym AUG Or thes experiment. ve and V3 


2 
Seeewnov Shown for clarity purposes as they have the same 
magnitude as Wms v, Poms lchcivyeerarecer and V3 Sebsecia c iby 
less than Yn Puuvewmewanape 15 identical. The units are 
maemeerary Since the basic equations are linear. However, 


W. 18 not plotted on the same scale as Vn: AiGUanci vies 


2 
em@ow COuble maxima with peaks at y = +1.2. These results 
are somewhat similar to those obtained by WSN without 
heating but their results showed the minimum between peaks 
as only slightly less than the maximum whereas in Figure 4, 
Peemminimum is less than one-half of the maximum. Aliso the 
Y1> Y23 and Yn fields are very close together but WSN 

showed Y4 To be predominate. The closeness of the yj fields 
roughly implies that the vertical mean streamfunction is 
Mearly zero so the fields are basically of a thermal charac- 
ter. Figure 5 contains the phase angles for Ya» Y2> Yn > 

and W., [Gmguumns experiment aS functions of y. v5) and Ym 

are very nearly coincident at the outer edges and va = Ym 
between y = te although only v5 is shown in Figure 5. V3 

is approximately 180 degrees out-of-phase from v5 and Ver 


which indicates the pressure trough at the surface is replaced 


rag 
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hig, saorowtun races for experiment 1. 
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- 4 6 8 10 


Pio 4 Simpdatudes sof Vins and w., for experiment 1. 


Ym and We are on different scales. 


EY, 








2; 
“180 


VW; 
w, 
Ge ses =e 
oe 
ol | | 
ia , 
“120 -60 0 60 120 180 


Jee ees Phases of V4 > V2> Wms and W for 


experiment 1 in degrees. Each dependent 
variable is expressed in the form 
Peos(kx+6), where P is the amplitude | 
and 6 is the phase angle. 





by a ridge at the upper level. This is consistent because 
the Ekman effect and heating are tied to the field at level 
3 which produces heating above surface low pressure systems. 
Uniake the results found by WSN for the adiabatic case, 
Pieure 5 shows the horizontal tilt of the disturbance to 
Bena Minimum at y = 0 and a maximum at y = +(1 to 4). The 
Bpeneral sense of the tilt is northeast but larger than that 
observed in easterly waves. Also the tilt of the ws, field 
memeLuates considerably. The small perturbations, particu- 
larly in the Wa field, suggest that there is not a clearly 
Gefined unstable mode in y. This does not agree with the 
meeutts of WSN which had no heating. It is interesting 
Gownove that the fields are very smooth but have large 
gradients except for W., Vie SMmOWS a radical behavior. 
Experiment 2 is the same as experiment 1 except the 
meemevetm 25 included. The results in general are the same. 
Figure 6 shows the growth rate profile. The magnitude of 
meemerowchn rates is almost the same although a small decrease 
in the maximum is noted (from 0.770 to 0.730) and the 
perturbations betwen ky, SO mnOmandul WO. have Deen somewhau 
smoothed. The maximum growth rate corresponds to a doubling 
time of 0.44 days and a wave length of 6280 km. Figure 7 
contains the amplitudes for experiment 2 as a function of 
Yee AS in experiment 1, only Ym and W. are shown with Vm 
and W, on different scales. Again there is a double maxima 
at y = +1.2 but the profile is asymmetric with the larger 


maximum in the shear zone to the north. This asymmetry is 
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Fig. 6 Growth rates for experiment 2. 
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Fig. 7 Amplitudes of Vn.and w, for experiment 2. 
Yen and W5 are on Geese reemr scales. 
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experiment 2 in degrees. Each dependent 
variable is expressed in the form 
Peos(kx-6), where P is the amplitude 

and 6 is the phase angle. 
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meecucameince the absolute vorticity gradient is not an 
mec, function of y. WSN found that the shift of symmetry was 
meeune south while Figure 7 shows a shift to the north. 


The phase relationships of yp Wm» and w, shown in 


ie? EE 2 
Figure 8 are very similar to those in experiment 1 and 
again only Y) is shown between y = ¢t2. 

Experiment 3 was run with the (4.1) wind profile and 
mmeruded CISK heating, surface friction, and horizontal 
Mmererion. ithe beta term was excluded. The horizontal 
Traction was added in an attempt to smooth out the vertical 
velocity (w,) pattern. Figure 9 shows the growth rate 
profile for experiment 3. The maximum growth rate has been 
reduced considerably (from 0.770 to 0.520) and the growth 
Moots off at KY G > 2.00. The maximum growth rate corresponds 
to a doubling time of 0.62 days and a wavelength of 6280 km. 
Note that there are no perturbations in the profile. Figure 
10 contains the amplitudes of Vn and W, Hore Ss experimenus 
Pema experiment 1, only Vm and W., are shown and are on 
different scales. Vy fomotiehely Varcer and 3 Slaaine by 
smaller than Yn: fhueercetaveresting To nove that the largest 
amplitudes are located near the boundaries for all fields 
whereas in the case without friction (experiment 1), the 
maximum amplitudes are in the central regions. Perhaps the 
large scale eddies favored by the horizontal friction are 
heavily damped by the horizontal shear between y = fil. 
Peeure 11 contains the phases of Vi> V2 and Ws Ym is not 


shown but it is everywhere equal to V5- Again Ym and V2 
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Fig. 9 Growth rates for experiment 3. 
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Piveawe so Oa v)> Yo, and W for experiment 


3 in degrees. Each dependent variable 
is expressed in the form Pecos(kx-6), 
where P is the amplitude and 6 is the 
phase angle. . 





are out-of-phase but all fields are smoother than in experi- 
ment 1. Where there are peaks in the amplitude fields, the 
phase angles of w, and 3 are the same which tends to 
suggest that warm air is rising and cold air is sinking. 
Again there is an eastward tilt in the J fields. The W,, 
field tilt is opposite to the w fields but it does correlate 
with the V2 field where the amplitude is a maximum. It 
appears the scale in y is about avi 

Experiment 4 is the same as experiment 3 but the beta 
term is included. The beta term has very little effect 
[eepe LO make Che solutions asymmetric. Figure le shows 
the growth rate profile which is the same as Figure 9 
Seecepy for the slight perturbation between the wave numbers 
0.40 and 0.90. The maximum growth rate corresponds to a 
doubling time of 0.62 days and a wavelength of 6280 km. 


ifmeeure 13 contains the amplitudes of Yin and w A Saat) 


5° 
experiment 3, vy and V3 are not shown. The only change noted 
mmoesliphnt increase in the maximum to the north and a 
decrease in the maximum to the south. Figure 14 contains 


the phases of wv and w, and again Yn is everywhere equal 


1? ¥3 
to v5 ° The W and Yo correlation is the same as in experiment 
S 

Experiment 5 was run with (4.2) wind profile and surface 
rection. Heating and horizontal friction were excluded. 


This profile is a mirror reflection of the UL 


is -U,. MimenorotmimMemicsamilar in general to the wind 


profile, that 


profile in the Intertropical Convergence Zone where there 
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is cyclonic shear in the lower region and anticyclonic 

shear in the upper region. Since the ITCZ is driven by 
heating and convergence, this is the type of mean wind 
pattern that is expected. Three cases were run although only 
one is shown in the table. Case 1, which is summarized in 
Table 1, had a grid mesh of 50 km, a y, value of 400 km and 

a time step of 1.0 hours. Case 2 used a grid mesh of 25 km, 
a i value of 200 km and a time step of 0.50 hours. Case 3 
used a grid mesh of 12.5 km, a Yo value of 100 km and a 

time step of 0.25 hours. Case 1 was run over the range of 
wave numbers from 0.10 to 1.20 and the maximum unstable 

mode was found at ky, = 0.65. Cases 2 and 3 were run only 

at wave number 0.65. Figure 15 shows the growth rate profile 
for experiment 5, case 1. The maximum growth rate at ky, = 0265 
was 0.283. The doubling time was 1.13 days with a wavelength 
of 3465 km. For case 2, the growth rate obtained was OAL ive 
with a doubling time of 0.85 days and a 1930 km wavelength. 
For case 3 a growth rate of 0.162 was found with a doubling 
time of 0.50 days and a 967 km wavelength. For cases @ and 
Ss rio 0.65 may not be the maximum growth rate. This shows 
that as the scale (y,) is decreased, the baroclinic process 
is suppressed. This is consistent Vaipmeparoeclinwe Iiscao— 
bility theory. It is noted that in case 3, the growth rate 
is the same as that obtained by WSN in their experiment l. 
Figure 16 contains the amplitudes of V4» Vn and w. for case 


f[eror clarity, !. is not shown but it is siimehtly Less 


6 
than Ya: These amplitudes are very similar to those of WSN 
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Fig. 12 Growth rates for experiment 4. 
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is expressed in the form Pcos(kx-6), 
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in their Figure 4A. The amplitudes for cases 2 and 3 are 
not shown but they are much the same as Figure 16. Figure 
17 contains the phases of V5» V2. Ym and W, as func Lons 


oe OY. Vn and w, are approximately 165 degrees out-of-phase. 


2 
mens implies a Daroclinic conversion of kinetic energy. 
meweyer, tne tilt 1S opposite to the shear which indicates 
MrTovrOnlececonversionm of kinetic energy. The tilt is 
westward with height north of the ITCZ and eastward with 
Mewene south of the ITCZ which is consistent with rising 
fereron ana Unsevtled weather that is often observed east 
mmmomencuriace trough axis in easterly waves north of the 
meen, (nese results are also consistent with those of 

WSN. 

Experiment 6 was run with the (4.2) wind profile, CISK 
meen ana surtace friction. The beta term and horizontal 
friction were excluded. Figure 18 shows the growth rates 
Pere Chis experiment with a maximum unstable mode at Kye =o le 
corresponding to a doubling time of 0.41 days and a wave- 
Heneth of 2285 km. The growth rate of 0.790 is similar to 
mae Of Experiment 1 and again the large value was due to 
the large value of n. The perturbations between wave numbers 
0.10 and 1.10 are even larger than in experiment 1 indicating 
meeeeconere iS no clearly defined unstable mode in y. The 
amplitudes have a pattern similar to Figure 4. The phases 


formebnas experiment tilted northeast as in the other cases. 
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VII. CONCLUSIONS 


Mmaniicomunestoan mumeri1calty computed barotropically un- 
stable disturbances are compared with tropospheric waves 
eat occur in the tropical easterlies. A two-level quasi- 
geostrophic model is used with an attached Ekman layer. The 
most unstable solutions are obtained numerically by the 
moeetal value procedure. Jwo basic shear zone profiles are 
used: (1) U, =U 


=-U. tanh WN and (2) U, =-U,=U, tanh WSIS 


5 3 
Calculations were made in which CISK heating and horizontal 
iMieeetlon were added. 

Certain features of the predicted wave structures do 
@errespond to observed tropical wave disturbances. The tilt 
of the disturbance in a sense opposite to the horizontal 
wind shear that was observed by Wallace and Chang (1969) is 
Meeoceie Jy) all the experiments. All experiments which contain 
heating have the warm structure which was observed by Chang 
ee al. (1970). The presence of rising motion to the east 
of the wave axis which was noted by Yanai and Nitta (1967) 
mmobserved in experiment 5. 

The growth rates obtained with heating are unrealistically 
large, but it was shown that these could be greatly reduced 
by choosing a smaller value of the parameter ne It was also 
noted that when heating is present there is no truely 
Meererrea Scale of instability. This has been observed by 


Weeme (1971). It was pointed out in the Introduction that, 


49 





while mature tropical wave disturbances may be driven by 
CISK heating, they may originate from barotropic instability 
Meaces 1970). If so, then phe @iecurbance should retain 
certain features (latitude structure, wavelength, etc.) 
of the barotropically unstable modes. This would provide 
a selection process which does not occur if CISK heating 
is present from the start. 

imemetidies of this thesis Should be continued but a 
more complicated model is needed. Proper parameterization 
of the heating term is necessary. A primitive equation 
model which includes a better treatment of the boundary 


layer should be used. 
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iemrepresented by U5 = U. = -U tanh Lae and the second is a mirror 


imme ction of the first and is given by U2 = -U,- The model is run 
for 142 days after which growth rates are determined as a function 
of the non-dimensional wave number Kyo: The associated phase 
angles and amplitudes for the wave numbers with maximum growth 
rates are given. The solutions with heating do not have a strongly 


Beeterred mode of maximum growth. 
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